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We review and extend the analogies between Gaussian pulse propagation and Gaussian beam 
diffraction. In addition to the well-known parallels between pulse dispersion in optical fiber and 
CW beam diffraction in free space, we review temporal lenses as a way to describe nonlinearities 
in the propagation equations, and then introduce further concepts that permit the description of 
pulse evolution in more complicated systems. These include the temporal equivalent of a spherical 
dielectric interface, which is used by way of example to derive design parameters used in a recent 
dispersion-mapped soliton transmission experiment. Our formalism offers a quick, concise and 
powerful approach to analyzing a variety of linear and nonlinear pulse propagation phenomena in 
optical fibers. 
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This paper introduces an ab-initio study of pulse prop- 
agation phenomena analogous to spatial CW diffraction 
behavior. We address both linear dispersive evolution 
as well the self-phase modulation effects of the nonlin- 
ear index of refraction |j| . The latter is responsible for 
much of the current interest in nonlinear optical com- 
munications, since pulse shapes such as solitons and 
dispersion-managed solitons display much more attrac- 
tive transmission properties than linear transmission for- 
mats (e.g. NRZ) || . 

Such nonlinear pulses are usually self-consistent eigen- 
solutions of a wave equation, which is the primary reason 
for their robustness to uncompensated spectral broaden- 
ing and resultant dissipation into the continuum. The 
conventional hyperbolic secant soliton is an exact solu- 
tion of the nonlinear Schrodinger equation || , and prop- 
agates indefinitely in a lossless medium without losing 
its shape. Lossless media can be realized in practice 
quite effectively by using lumped amplification stages, 
and erbium-doped fiber amplifiers offer excellent charac- 
teristics in this regard. 

Breathers, sometimes called dispersion-managed soli- 
tons H, are also self-consistent 'eigen solutions' of 
the wave equation that propagate with periodic pulse 
width, chirp etc. While not strictly unchanging in shape, 
breathers evolve back to their initial configuration, essen- 
tially traversing a closed, non-degenerate orbit in phase 
space [||. Unlike pulse shapes designed for linear trans- 
mission channels, these pulses do not require periodic 
dispersion compensation along the transmission channel, 
and so offer an attractive alternative to the strong control 
requirements of the nonlinear Schrodinger soliton. 

Characterizing the solutions of the nonlinear wave 
equation is often simplest via direct numerical simula- 
tion, and this has been particularly true for dispersion 



mapped solitons 0. In order to understand, capture and 
then predict and utilize the essential physics that guides 
this behavior, a more conceptually accessible framework 
is sometimes preferable, such as the variational approach 
with a pulse shape Ansatz ||. The pulse shape is de- 
scribed as a dynamical system; we write the Hamilto- 
nian based on the action principle and seek solutions to 
the Euler-Lagrange equations of motion |], [l(J . This ap- 
proach is not always applicable, however, especially when 
the Ansatz is incapable of capturing some essential phys- 
ical behavior. Also, it is somewhat more of an analytical 
tool for probing the dynamics of systems that we already 
know something about, or can predict at least partially, 
and it may be convenient to have other approaches that 
can offer quick insight into constructive aspects of non- 
linear propagation, so that different geometries can be 
analyzed and compared quickly and easily. 

The parallels between dispersive pulse propagation in 
optical fibers and paraxial CW Gaussian beam diffraction 
in free space have been identified for some time jll], [l^, 
[D|. More recently, the analogies have been extended to 
include temporal lenses as a way to translate the imaging 
properties of spatial lenses into the temporal domain p4[ . 
In this way, pulse correlation and convolution devices 
may also be constructed 15 1. Still more recently, it was 



shown that temporal lenses can characterize nonlinear 
effects in the wave equation, leading, for example to the 
formation of a class of steady-state repeating pulses jl6) . 
We believe that this is perhaps the most potentially use- 
ful of the space-time analogies: in this paper, we further 
extend the use this formalism to describe still more pow- 
erful applications such as Gaussian pulse propagation in 
optical fiber systems, including dispersion mapped sys- 
tems, including the effects of the nonlinear index of re- 
fraction. 
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We first outline the basic physics that motivates this 
discussion and sets the context for further development. 



2 



I. SPACE-TIME ANALOGY OF BEAM 
DIFFRACTION AND PULSE PROPAGATION 

A. CW Gaussian beam diffraction 

The Fresnel-Kirchoff diffraction integral is a well- 
founded approach to electromagnetic propagation prob- 
lems, and several textbooks cover the topic from a variety 
of approaches 0, [l^, . We will briefly review only as 
much as necessary to establish our argument, limiting 
our argument to diffraction in f +1 (x,z) dimensions. 

An electromagnetic field of radian frequency u> and 
scalar complex amplitude u(x, z) can be represented 



E(x,z,t) = u(x,z) exp(zwt) 
where u(x) obeys the wave equation, 



(1) 



V 2 u + k 2 u = 0, fc 2 = u V=fy) ■ (2) 

This equation admits plane wave solutions of the form 
exp(±ifcz) representing propagation along =Fz respec- 
tively, and indeed, an arbitrary superposition of plane 
waves, each with the same wavelength, propagating along 
all possible directions, 



(x,z) = / u (k x ) exp[i(k x x) -i^k 2 - k 2 z\ dk x (3) 



where uq is the Fourier transform of the input field 
u o (x,0). 

We consider optical beams whose plane wave compo- 
nents propagate at small angles to the z axis (paraxial 
approximation), so that we can expand the square root 
in teh in a Taylor series and keep the first two terms, 



E(x, z) 



kz J u (k x )exp (i^fcz^J exp(ik x x) dk 



(4) 



-iojt—ikz 



'£-Ju {x')exv 



ik(x-x') 2 
2z 



dx' 



where the term in parentheses defines u(z,t), the field 
envelope, 



u(z, t) 



ik 
2ttz 



uq(x') exp 



ik 

Yz 



>\2 



(x-af) 



The propagation of continuous-wave (CW) Gaussian 
beams in free space and rotationally-symmetric quadratic 
graded-index media is conveniently described by assum- 
ing that the envelope has the form fl7f 



exp 



P(z) + 



2q(z) 



(6) 



where, we find by substitution into the wave equation (|2j) 
that dP/dz = —i/q(z) for such media. The g-parameter 
describes the Gaussian beam completely, 

q(z) R{z) Trnw 2 {z)' { ' 



In the above definition, R{z) describes the radius of cur- 
vature of the beam, and w(z) the beam spot size. 

The usefulness of the q-parameter lies in the bilinear 
transformation (ABCD law) that characterizes how this 
parameter evolves with propagation. For an optical sys- 
tem described by a real (or complex) ABCD matrix, the 
output q parameter is given by 



q 



Aqj + B 
C qi +D' 



(8) 



Separating the real and imaginary parts of q D enables 
us to calculate the radius of curvature and spot size 
of the Gaussian beam at the output of the optical sys- 
tem. Many practically important optical systems and 
their corresponding phenomena can be described by sim- 
ple ABCD matrices, such as propagation in a uniform 
medium, focusing via a thin lens, beam transformation 
at a dielectric interface, propagation through a curved di- 
electric interface and thick lens, propagation in a medium 
with a quadratic index variation etc. Table 2-1] 



B. Gaussian pulse propagation 

Consider a single mode in an optical fiber, usually the 
lowest-order fundamental mode, excited at z = 0, and 
with an assumed temporal envelope of the Gaussian form, 



E(z = 0, t) = Re [exp(-crf 2 + iw t)] 
and write as a Fourier transform integral, 



E(0,t) = Re 



exp(iu>ot) / uq(Q) exp(iilt) dtt 



(9) 



(10) 



where Uq is the Fourier transform of the Gaussian enve- 
lope Mo = exp (—at 2 ). 

As in the spatial case, we can be think of this as a su- 
perposition of time-harmonic fields, each with frequency 
(ljq + Q,) and amplitude ii (Q)dfl. These waves will expe- 
rience a phase delay when propagating a distance z; we 
multiply each frequency component by its propagation 
delay factor exp[— ifl(u>Q + Sl)z] so that 



dx' (5) E{z,t) = J uo(Jl) 



fill 



Expanding (3(u>q + f2) in a Taylor series about the op- 
tical frequency u)q, 



0{u o +a) = p(co ) 



d£_ 

du> 



9. 



1 (P(3_ 
2 



we can write 

E(z,t) — exp[i(w £ — (3qz)] x 
uo(O) exp <, i 



Vlt - (3'Qz - -(3"n 2 z 
2 



(12) 
(13) 

• dn 
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where the integral defines the field envelope u(z,t), so 
that 



E(z,t) = u{z,t) exp[i(ujQt — Pqz)]. 



(14) 



The differential equation satisifed by u is, to second order 
of derivatives of || , 



du ,du 1 d 2 u 



(15) 



The solution (|13| ) can be written using the inverse Fourier 
transform relationship, 



u(z, t) 



u Q (Q) exp (-^f3"fl 2 z - i0'Qz^ 
x exp(iSli) dfl 



1 



i2n0"z 

x / uo(t') exp 



2(3"z 



— (T - t'y 



dt' (16) 



where T = t — 0'z = t — z/v g is the time coordinate in 
the frame of reference co-moving with the pulse envelope 
at the group velocity v g — 1//3'. Dispersion of the group- 
velocity (GVD) is represented by 3". 

The formal similarity between (^J) and (|l6| ) is the prin- 
cipal motivation for this analysis. We can write down 
a set of space-time translation rules (see Table |) to ap- 
ply results from spatial diffraction to temporal dispersion 
and vice versa. One family of results that can be derived 
from this space-time analogy correspond to spatial imag- 
ing e.g. the 2-f and 4-/ optical systems. These can be 
applied to pulse compression or expansion experiments 
etc @. 

But we will see in later sections that many linear and 
nonlinear pulse propagation systems can be described by 
cascading simple ABCD matrices, and this can result in 
substantially simpler calculations and more direct phys- 
ical understanding of the physical processes involved in 
nonlinear pulse propagation. We will first need to develop 
some additional facility in characterizing optical systems 
associated with the pulse propagation equations. 

The spatial g-parameter has a temporal equivalent q t 
in accordance with the space-time translation rules of 
Table |, defined by 



1 



q t {z) R t (z) 



+ i 



28" 



(17) 



where r(z) represents the pulse width (scaled in the T 
frame by V2) and Rt(z) its chirp. A Gaussian pulse 
in linearly dispersive fibers is then represented by the 
envelope fllq] 



u(z, t) 



u 



to 



r(z) 



(18) 



exp 



i tan \- i 



3" 



Co 2B"R t {z) \B"\ t 2 (z) 



where the pulse width and chirp satisfy evolution equa- 
tions in linear dispersive fibers exactly analogous to their 
spatial counterparts, beam spot size and radius of curva- 
ture, in free space Era 



(19) 



t\z) =7^(1 + ^), 
R t {z) = 2(1 + % 



with Co =Tq/2|/3"| defining the dispersion length 



II. COMPONENTS OF THE ABCD 
FORMALISM FOR GAUSSIAN PULSE 
PROPAGATION 

As a simple example of the application of the above 
translation rules, we consider the propagation of a Gaus- 
sian input pulse with envelope 



£7(0, T) = exp 



T 2 
2T 2 



(20) 



The transmission medium comprises of two concatenated 
sections of fiber with lengths z\ and z^ and with GVD 
coefficients 3" and 3'^ respectively. We ignore any non- 
linear effects in this simple problem, and assume that the 
medium is lossless. What are the pulse characteristics at 
the output of the second medium i.e. what is the pulse 
width at z = z± + Z2? 

One way of solving this problem is by recourse to the 
wave equation ( |l6[ ) solution by the Fourier transform 
technique. We have, 



U{z2,w) = U(z u u)exp ( u -3'iz x u 2 



U(z ,u>)exp 



{B>{z 1 + 8'iz 2 )uj< 



.(21) 



Taking the inverse Fourier transform, 



U( Zl +z 2 ,T) = 

i r 00 ~ 

5- / U(0,w) 



- {B"zi + 82 z 2 ) lo 2 - iuT 



T 2 



du> 



T 2 - i{B'{ Zl + 8^2) 
T 2 

x exp 



2(T 2 -i(0'{z 1 + B 2 !z 2 ) 



(22) 



from which we see that the ratio of the output to input 
pulse width, therefore, is 



Tx 
To 



flZi + B 2 Z2 



T 2 
J o 



(23) 
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TABLE I: Space-time translation rules 



spatial frequency (Fourier variable) 


k x 


Q 


frequency (Fourier variable) 


transverse distance 


X 


t--*- 


time (in moving reference frame) 


propagation distance 


z 


z 


propagation distance 


wavevector (inverse) 


k- 1 


-P" 


GVD coefficient (negative) 



We will now verify (|23|) using the ABCD matrix ap- 
proach. The system is described very simply by the prod- 
uct of three matrices, 



M = 




so that 



92 = 



Agx+B_ _ , % 

Cqi+D 



Z 2 - 



Using the shorthand notation 



R 2 = R{z 1 + z 2 ) 1 R x ee i?(0), 

t 2 = t(zi+z 2 ), TlEEr(O). 



(24) 



(25) 



(26) 



we have 
i 

B.2 



ft 



1 

Rx 



; 2/3J' 



(±)+(f) *(*)'♦(*) 



M Z1 



Z-2 



(27) 



The real and imaginary parts of both sides of the above 
equation have to be equal, leading to a pair of simulta- 
neous equations. For an unchirpcd input pulse, R\ = 
so that equality of the imaginary parts leads to 



1 

~R~> 



2f% 



\T 2 Tl 



Subsituting this expression into the equation of equality 
of the real parts of (|27j ) and some algebraic manipulation 
leads to 



t(zi + Z 2 ) _ 7£ 

r(0) ~ n 



ftZ! + fj' 2 'z 2 
rf/2 



(28) 



which is the same as (^3|), since r 

In the above calculation, we have carried out some 
algebraic simplifications by hand in order to show that 
the result obtained by the ABCD matrix approach is the 
same as that obtained by the Fourier transform approach. 
Nevertheless, the former is computationally much sim- 
pler, and separating the real and imaginary parts of (|27 



as part of a numerical algorithm can be carried out with- 
out the notational complexity of, for example, rational- 
izing the denominator. 

While second-order dispersion is conveniently repre- 
sented by the ABCD matrix approach, there are prob- 
lems with extending the analysis to higher orders of dis- 
persion. The slowly-varying envelope equation analogous 
to ( JTsj ) including the effects of third-order dispersion 



au _ i a*u 



dz 



Q T 2 



6 P dT 3 



(29) 



or its solution in terms of the Fourier transformed vari- 
ables, 



U(z,u>) 



U(0,uj)exp ( ~/3"zu 2 + qP"'zu 3 



(30) 



does not have an equivalent in the CW spatial diffraction 
context. 

To see this, consider the next term in the Taylor ex- 
pansion of yjk 2 — fc 2 , in (||), which leads to an expression 
of the form 



u 2 {x) = e 



—ikz 



(31) 



r / £2 k 4 \ 

/ u^k^exp ii-^z + i-^zj exp(ik x x) dk x . 

Using the space-time translation rules, we find that the 
above expression contains a description of second and 
fourth-oidev dispersion, not third-order dispersion. 

This is obviously a general characteristic of the above 
Taylor expansion; all odd-order dispersion terms have no 
spatial paraxial diffraction equivalent in the ABCD ma- 
trix content. Recall that the effect of j3' is accounted for 
by transforming to a moving reference frame T = t — [3 1 z. 

For completeness, we derive the translation rule for any 
even-order dispersion in terms of the equivalent term in 
CW diffractive optics. A little algebra will show that the 
generalization of (E30|) yields 



U(z,u) = U(0,u) x 



+ --- + tAt/3 (2 '' ) ^ 



(2r) 



(32) 



(33) 



and, correspondingly, for the diffraction of a Gaussian 
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beam, 



U(x,k x ) = U(0,k x )exp 



l 2k Z + l W Z ' 



+l- 



r\ J -- L V2 



„2r 



k 2r-l' 



(34) 



Therefore, the translation rule for 2r-order temporal dis- 
persion is given by 



2r) 



(2r) 



-! 2' 



rs " \ 2 



n(2/-i) 

i=i 



Note that this corrects the statement in |l 

The slowly varying envelope equations corre- 
sponding to modulated plane waves in disper- 
sive media have the same form as the parax- 
ial equations describing the propagation of 
monochromatic waves of finite spatial extent 
(diffraction). 

We append that this correspondence holds for all even 
orders of dispersion, and of course, for (3' as well, by 
transforming to a moving reference frame. 

Our ABCD formalism would be of limited interest if 
the only phenomena it could capture were that of disper- 
sive propgation. But, as mentioned in an earlier section, 
the development of the time-lens formalism lets us de- 
scribe nonlinear mechanisms as well. 

By analogy to spatial lenses which are characterized by 
a lens factor exp(ifcr 2 /2/) which multiplies an incoming 
optical beam, we define a temporal lens as a device that 
multiplies the pulse envelope by a factor (T|, [HJ 



Lens Factor = exp 



if- 



exp 



-ibt 2 ] 



(36) 



The ABCD matrix representing a temporal lens has 
the same form as that of a spatial lens, 



M 



1 




(37) 



where ft represents the temporal "focal length" . 

A comparison of spatial and temporal lensing is shown 
in Figure |l| In the spatial case, the lens compensates for 
the spreading of the beam waist, and "flips" the phase 
fronts to convert a diverging beam into a converging one. 
Similarly, a temporal lens reverses the sign of the chirp, so 
that further propagation in a j3" < dispersive fiber will 
compensate for the chirp (phase modulation) caused this 
far. This is also an interesting and physically illuminating 
approach to discussing the physics of the formation of 
solitons |l7|, Chapter 19]. 
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FIG. 1: (a) Spatial lens (b) Temporal lens 



One possible implementation, as proposed in |T(| , is to 
achieve temporal lensing by self-phase modulation during 
the passage of the pulse through a section of nonlinear 
fiber (/?" » 0, n 2 > 0). For short distances, z < ttt^/\P"\ 
and when /3"/ r o ^ (27ro2/A)I p for peak intensity I p , a 
pulse with input electric field envelope u(0,T) emerges 
from a length z of nonlinear fiber with phase modulation 



u(z, T) = u(0, T) exp 



-i— |«(0,T)| 



2ctj 



(38) 



where r\ = yjife defines the impedance of free space. If 
we write the pulse intensity as 



— =/ p exp 



-2,1 

TO 



(39) 



and keep the first two terms in the Taylor expansion of 
the exponential in (pq) . 



u(z, T) = u(0, T) exp 



. 2u; n 2 IpZ 



(40) 



modulo a phase term linear in z that is independent of T. 
The effect of propagation through length L of nonlinear 
fiber is to impart a quadratic chirp to the pulse, which we 
represent by the multiplicative term exp(— ibt 2 ) so that 



1 



2ui n 2 IpL 



W'ft 



(41) 
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Another method of obtaining time lensing is based on 
the principle of electro-optic modulation . An electro- 
optic phase modulator driven by a sinusoidal bias voltage 
of angular frequency uj m results in a phase modulation 
that is approximately quadratic under either extremum 
of the sinusoid. The phase shift can be written as 



exp[i(f>(t)] = exp 



-iK 1 



2 +2 \ 1 



(42) 



where K is the modulation index |l7], §9.4]. In this case, 

2 



Kcul 



(43) 



We have described our temporal lens by a section of 
nonlinear fiber of j3" « 0, analogous to a spatial thin lens, 
which is assumed to have no thickness. Just as practical 
lenses do have some thickness, practical fibers have non- 
zero (3" . For those situations in which this cannot be 
ignored, or may even be utilized constructively, we derive 
the corresponding equivalent of a spatial "thick lens" . 

Our first step is to characterize the temporal equivalent 
of a curved dielectric interface: a spatial lens comprises 
of two such interfaces separated by a length of material 
of enhanced refractive index. At a planar dielectric in- 
terface between two media of refractive indices n-y and 
n,2 , a Gaussian beam undergoes a change in the radius of 
curvature, but is unchanged in beam width, 



i?2 — —Ri, 
m 



W\ = u> 2 - 



(44) 



By analogy, a chirped Gaussian pulse at the interface be- 
tween two fibers of GVD coefficients f}'{ and 0% trans- 
forms to a different chirp, but with unchanged pulse 
width, 



f 



1 



/?2 R2 A" Ri 



n = r 2 . 



(45) 



Of course, the pulse width evolves differently in the two 
sections of fiber, 



'CM 



z 2 



7=1,2 



(46) 



where £oi is the dispersion length in fiber i. 

The ABCD matrix for a (spatial) spherical dielectric 
interface and its temporal translation are 



M 



1 







rt 2 — «i rti 
n 2 R ri2 



1 \ 

Ri Pi' ) 



What does this represent? We use the ABCD bilinear 
transformation, 



92 = <2i J 
which implies that 



1 

?2 



1 

R~2 



1 

Ri 



02 \ 02 

™ ) Qi + — 



(48) 



Ri \ % 



1 2/3" \ /3'J , Jn , 



After some algebraic manipulation, we can write the 
above as 



1 

Ri 



1 

R~2 



l 



1 



Ri R\ 



(50) 



showing explicitly how the chirp transforms at this inter- 
face. 

The ABCD matrix for a temporal lens of "thickness" 
d is written as the product of three ABCD matrices rep- 
resenting, when read from right to left, a transition from 
the input fiber to the fiber that defines the thin temporal 
lens, propagation in the second fiber, and a transition 
back to the input fiber, 



M = 




1 

zKM. 

Ri 



and, multiplying the matrices together, we get a single 
ABCD matrix which defines the output q t parameter via 
the usual bilinear transformation (Aq t + B)/(Cqt + D), 



M = 



-Ra) 



(f -0(> 



1 + ^L(l_4 



02 . 



d 

-R 2 



The temporal focal length f t is analogous to the spatial focal length and is given by — A/C, 



ft = 



Ri \ 0'i 



l 

R~i 



1 



1-^]- 

0'i 



Ri(-R 2 ) 



0i 



Si 



1-% 



(51) 



(52) 



The temporal focal length defines the time from the out- 



put plane at which an initially unchirped pulse becomes 
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unchirped again. 

We can write the above in slightly simpler notation, 
for the specific case R± = — R2 — R, and let k — d/R, 
A/3" = ft - ft, 



1 — AC 



A/3" 



k A/3" A/3" 



i? 
~2 



(53) 



where the term in parentheses represents an enhancement 
factor over the "thin lens" formula. 

For k <C 1, we can simplify the above expression keep- 
ing terms of O(k), 



1 



1 



k A/3" 
2 



1 



A/3"\ 2A/3" 



k A/3" ' 

2 7T; w 



2A/3" 



(54) 



The above relation confirms our physical intuition that if 
ft — ft — A/3" < 0, then we have reduced ft, the distance 
to the point of zero chirp from the output plane, for an 
initially unchirped input pulse. 

We now have the tools we need to analyze a reasonably 
complicated practical problem: designing the length of a 
dispersion map so as to get self-consistent eigen-pulses 
with periodic pulse width and chirp. 



III. DISPERSION-MANAGED SOLITON 
TRANSMISSION EXPERIMENT 

It has been recently found that a stable, self-consistent 
pulse solution exists in a dispersion-managed fiber trans- 
mission system ||. While these are not solitons in 
the strict mathematical sense, they have been called 
dispersion-managed solitons, or perhaps more appropri- 
ately, breathers. They demonstrate periodic behaviour: 
the pulse width and chirp of Gaussian breathers, for 
instance, are periodic functions of the propagation dis- 
tance. Breathers share a property in common with soli- 
tons in that they can propagate indefinitely without 
losing shape; even though the pulse shape undergoes 
changes within a disperion map period, the pulse does 
not disperse away to infinity, or tend to self-focus to a 
point either of which invalidate the applicability of the 
nonlinear Schrodinger equation after a certain distance. 

A dispersion-mapped (DM) soliton is closer to a Gaus- 
sian shape than the hyperbolic secant of the nonlinear 



Schrodinger equation ]2Q] , and it is interesting to ask 
whether our analysis is capable of capturing the essen- 
tial aspects of its evolution along a dispersion-mapped 
transmission channel. 

We consider, as our example, the paper by Mu et 
al. |HJ who have simulated DM soliton dynamics in a 
recirculating fiber loop. Their dispersion map consists of 
100 km of dispersion shifted fiber (SMF-LS) with nor- 
mal dispersion D\ equal to -1.10 ps/nm-km at 1551 nm, 
followed by an "approximately 7-km span" of standard 
single-mode fiber (SMF-28) with an anomalous disper- 
sion D2 equal to 16.6 ps/nm-km at 1551 nm. The re- 
sults of the paper indicate that Gaussian shaped pulses 
of pulse duration 5.67 ps and peak power 9 dBm were 
used. We will derive the result that, for these parame- 
ters and given the length of SMF-LS fiber, the length of 
SMF-28 fiber that needs to be used is indeed "approx- 
imately 7-km" . In other words, we will show that this 
given dispersion map can support lowest-order chirped 
Gaussian self-consistent solutions, i.e. breathers. 

The dispersion map, shown schematically in Figure |^, 
consists of three fiber segments: a length z\/2 equal to 
50 km of SMF-LS fiber, followed by a length z 2 of SMF- 
28 fiber, whose numerical value is to be determined, and 
then the remainder Zi/2 of SMF-LS fiber. Each segment 
of fiber has nonlinear characteristics, which we model via 
a time lens situated, for simplicity at the individual mid- 
points of the respective segments. Consequently, each 
segment is described by the cascaded product of three 
ABCD matrices, with two additional matrices represent- 
ing the transitions between fibers of different 0". For 
simplicity, we will assume that the nonlinear properties 
of the fibers are identical. 

The overall ABCD matrix for the system can be writ- 
ten down quite easily, 



M 




which, after some algebra, can be written as an ABCD 
matrix with the following elements, 



8 



.4 



D = — -i 



Zl_ 

4/ 



V 







v 4/ 
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_*\ 


4/ 


2/y 



2 - 



«2 

2/ 



f2 

2/ 



V 4/ 



£l 
4/ 



1 



2- 



4/ 



2-^ 
4/ 



4/ 



2/ 



4/ 



2/ V 2/ 



VJ V 4/ 



(56) 



1 ~ Tf) ( " 
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i - 
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1-^ 
4/ 



2- 



2-^ 
2/ 



£2 

2 



£2 

2/ 



1-^ +^1-^ 



4/ 



4/ 



£l 
4 



1 - -4: 



2-^ 
4/ 

2/ 



2/ " 



4/ 



£2 
2/ 



£2 
2/ 



1 - -r 



4/, 



(58) 
(59) 



The algebraic complexity of writing out the expressions 
explicitly should not mask the simplicity of multiplying 
two-by-two matrices, usually numerically. Note that the 
expression (|57]) for D is algebraically identical to that for 
A ([36]) , and it may be verified that AD — BC = 1 . 

The g-parameter (we have dropped the t subscript in 
this section for notational elegance) evolves according to 
the bilinear transformation law, and we require that the 
pulse repeat itself after propagation through one such 
ABCD matrix, 



1 A + B/g 
q~ C + D/q 

which has the solution 



(60) 



1 D- A 



1- 



q 



2B 



±i- 



D + A 



B 



(61) 



Since D — Am our above analysis, we already see that 
q is purely imaginary at z — i.e. the pulse has zero chirp 
at the midplanes, as we would expect a breather to have. 

At this stage, we can substitute numerical values for 
the various parameters (except zi , which is what we seek) 
into the expressions for the A, B, C and D elements (|5^ - 
p7| ) and solve ( |6l| ) numerically for z%. While this is not 
difficult, and already yields a quick solution to the prob- 
lem at hand, we can get further insight via a well-justified 
simplification as follows. 

The q parameter at the midplanes, where it is purely 
imaginary, is given by 



1 



2|fl'| 



(62) 



where j3" = 1.40 x 
r = 5.67 x 10 



10 



-27 „2 



s /m and input pulse width 



12 s. Consequently, for such pulses, 1/q 
0, and since A = D, this implies that A = 1 in (|6l|). 



With the notational substitutions 



£2_ _ fl_ all _ Pl_ 

X 2/' V 4/' Pr % 



(63) 



we get the necessary condition 

[(1 - yf - (2 - y)y] (1 - x) - fly (2 - y)(l - y) 

-^(l-y)(2-x)x = l. 



The solution of this equation is given by 



,,2/(2-2,) 



x — (3 

i -y 

or, in terms of the initial variables, 



(64) 



(65) 



Z-2 



2 4/ 



4/ 



(66) 



which is the necessary condition in order to have a stable 
self-consistent Gaussian eigen-pulse (breather) solution 
to the dispersion-map problem. 

All that remains is for us to interpret the variables in 
terms of the original problem and numerically evaluate 
this expression to get the desired length z 2 of SMF-28 
fiber in this dispersion map. The various numerical val- 
ues are as follows: 



TO 
Zl 



= 1.40 x 10~ 27 s 2 /m, 

= -2.12 x 10~ 26 s 2 /m, 

= 5.67 x 10~ 12 s, 

= 10 5 m 



9 



SMF-LS 



SMF-28 



SMF-LS 



IV. HERMITE-GAUSSIAN BASIS 




zy 2 



(a) 




Z1/2 



V / 



time lens 



GVD transition propagation 



(b) 



FIG. 2: (a) Analytical schematic of dispersion map from |2l]] 
and (b) its representation to express in terms of ABCD matrix 
elements. 



Our ABCD matrix formalism for pulse propagation ap- 
plies to chirped Gaussian pulses. To analyze more com- 
plicated shapes, we can expand the given pulse shape in a 
basis of chirped Hermite-Gaussian functions, which form 
a complete orthonormal basis (l^, ^0). The Hermite- 
Gaussian function (we consider only unchirped Gaussians 
here for simplicity) of order n is denned as the product 
of the Hermite polynomial of order n with a Gaussian 
function, 



if> n (t) = H n {t)exp(-t 2 /2), 



(68) 



where, for example, 

H (t) = l, Fi(*) = 2t, H 2 (t) = 4t 2 -2. (69) 

We can expand an arbitrary input amplitude itrj(i) in 
this basis, analogous to expanding a field in terms of 
plane wave components, as in solution techniques of the 
standard parabolic diffraction equation by means of the 
Fourier transform, 



M oW = ^c„H„(t)exp(-t 2 /2) 



(70) 



Given the nature of the problem, we realize our time lens 
with the nonlinear fiber as described earlier (^l|) , so that 



1 

7 



NL 



2n^\ n 2 I p L NL 
T 



(67) 



and take (3% L = f%, I p = 3.62 x 10 6 W/M 2 so that with 
fiber core area A c g — 47 fim, we get P = 8 mW = 
9 dBm. Also, we take Lnl — z 2 consistent with our 
choice of (3'^ L . 

The numerical solution (of the quadratic equation) for 
z 2 is equal to 7.00 km which is indeed the value "approx- 
imately 7 km" stated in the paper . In spite of appar- 
ent exact agreement, we should be careful to appreciate 
that this analysis is a characterization of only the most 
important processes in this experiment. Possible sources 
for approximation include the fact that a DM soliton is 
only approximately Gaussian, and that we have repre- 
sented the combined dispersive and nonlinear properties 
of the fiber segments by a single temporal lens. A better 
approximation may be to include several temporal lenses 
for each segment of fiber; this would make the algebraic 
expressions in this paper quite cumbersome to write down 
explicitly, but the numerical computation would not be 
much more difficult, since the matrices are only two- by- 
two in size, and comprise of purely real elements. The 
experimental configuration of pl| also includes several 
other elements which can affect the pulse shape, such as 
filters, fiber amplifiers and polarization controllers. 



where because of orthogonality of the Hermite-Gaussians, 
the expansion coefficients are given by 



1 



/7r2"n! 



u (t)H n (t)eM-t 2 /2)- 



(71) 



The propagation equation (|T^) defines the output pulse 
shape as the convolution of the input shape with a Gaus- 
sian kernel. Hermite-Gaussians, when convolved with a 
Gaussian, yield the product of a Hermite polynomial and 
a Gaussian p8[, 



dt ip n (to) exp --(t-t ) 2 



at 2 - 



2tt fa- I 



a+ 1 V a+ 1 



exp 



2(a 2 - 1) 



1pn 



vV-1 



(72) 



Taking as input the n-th Hermite-Gaussian mode uq — 
i/} n (t) (which has width To = v2), we evaluate the am- 
plitude of this mode after propagation through distance 

z, 



u n (z,T) = 



1 



x exp 



1 + ij3"z 
i T 2 



1 



0»z 



2(3" z 1 + 



WW 



(i"z 



n/2 



T 



(73) 



which can be seen to agree with (|l9|). 

A Hermite-Gaussian therefore maintains its shape dur- 
ing propagation, but adds a chirp (which is the same for 



10 



all modes) and a scaling of the width according to (19). 
Power conservation implies that the amplitude corre- 
spondingly scales down. The only term that is dependent 
on the order of the Hermite-Gaussian is a phase term; 
higher-order modes have greater phase advances, since 
their spectral content is higher. The important observa- 
tion is that the orthogonality of the Hermite-Gaussian 
expansion is preserved, and so this expansion may be 
used to predict the pulse shape obtained by propagating 
an input pulse. Our formalism remains valid as long as 
the differential equation describing the propagation of a 
particular order Hermite-Gaussian is of the form (15), 
i.e. the slowly-varying envelope approximation is valid. 
Therefore, we can expect that the lower-order expansions 
are usually valid; the results of applying our analysis to 
higher-order expansion terms generate the residual field 
corrections to the lower order results 11231 ■ 



V. CONCLUSIONS 

We have developed a 2 x 2 ABCD matrix formalism 
for describing pulse propagation in media described by 
Maxwell's equations, accounting for dispersion, nonlin- 
ear and gain/loss mechanisms. The method is analo- 
gous to techniques used in CW beam diffraction analy- 
sis, and correspondingly similar phenomena can be pre- 
dicted, such as chirp transformation, focusing, periodic 
pulse width expansion and narrowing etc. The spatial q 
parameter has a time equivalent qr in accordance with 
the given space-time translation rules. The real and 
imaginary parts of q^ 1 represent the chirp and the width 
of the pulse as a function of propagation distance z. 

The propagation of various input pulse shapes can be 
described by expanding the given pulse in a basis of 
Hermite-Gaussian functions; the ABCD formalism ap- 
plies to each Gaussian wave function separately. Propa- 
gation through a complicated system of optical elements 



is simple to calculate in terms of ABCD matrices: the 
resultant matrix is the cascaded product of the ABCD 
matrices of each of the individual elements with the ap- 
propriate ordering. The overall qx parameter is given by 
a bilinear transformation in terms of the ABCD elements 
of the overall product matrix, exactly analogous to the 
spatial case. 

We have formulated ABCD matrices for pulse propa- 
gation in dispersive fibers, and for temporal lenses which 
can characterize self-phase modulation phenomena. A 
spatial dielectric interface translates to an interface be- 
tween fiber segments of dissimilar GVD coefficient j3" . 
The temporal equivalent of a curved dielectric interface 
is useful for characterizing the transition between such 
dissimilar 0" fibers with the added presence of fiber non- 
linearities arising from the nonlinear index of refraction 
7i2- We have used these tools to characterize a reason- 
ably complicated real-life system: calculation of the dis- 
persion map for self-consistent stable propagation of a 
dispersion-managed soliton. 

We believe this method of analysis forms a useful 
complement to conventional pulse propagation methods, 
such as the split-step Fourier transform numerical proce- 
dures which are substantially more computationally 
intensive. The ABCD approach is useful for clarifying 
the important dispersive and nonlinear focusing effects 
in dispersion-mapped nonlinear fiber segments. Together 
with the variational approach S , based on modeling the 
pulse as a dynamical system characterized by a Hamilto- 
nian functional [^ , the g-parameter offers an insight into 
pulse evolution from a theoretical standpoint. 
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